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Abstract 

This paper is concerned with the asymptotic behavior of the solution to the following 
semilinear parabolic equation 

ut - At* + /(«) = 0, (i,t)e!lxR + , 

subject to the dynamical boundary condition 

d v u + fiu + ut = 0, (x,t)eTxR + , 

and the initial condition 

u\t=o = Uo(x), x G fl, 

where f2 C M. n (n £ N + ) is a bounded domain with smooth boundary T, v is the out- 
ward normal direction to the boundary and \x £ {0, 1}. / is analytic with respect to 
unknown function u. Our main goal is to prove the convergence of a global solution 
to an equilibrium as time goes to infinity by means of a suitable Lojasiewicz-Simon 
type inequality. 

Keywords: Semilinear parabolic equation, dynamical boundary condition, Lojasiewicz- 
Simon inequality. 

Introduction 

this paper we consider the following semilinear parabolic equation 

u t - Au + f(u) = 0, (x,t)ettxR + , (1.1) 



subject to the dynamical boundary condition 



d v u + fiu + ut — 0, 



(x,t) e T x R+ 



(1.2) 



and the initial condition 



U\ t =0 = U (x) 



iGSI, 



(1.3) 



In the above, Q C R n (n G N + ) is a bounded domain with smooth boundary T. v is the 
outward normal direction to the boundary and /i G {0, 1}. 

Parabolic equation and systems with dynamical boundary conditions have been ex- 
tensively studied in the literature (see for instance, [3lffll71-[9|fT2 | [T5 | [T9] etc.). In particular, 
local existence and uniqueness of solution to general quasilinear parabolic equation (sys- 
tems) with dynamical boundary condition has been established in a series of papers by 
Escher [7H9] (see also |15| for a semigroup approach in the ifp(fi)-setting, and jl] for the 
solvability result in a weighted Holder space). Moreover, in [8j|9], the author showed that 
the solution defines a local semiflow on certain Bessel potential spaces including the stan- 
dard Sobolev space W 1,p . Based on the approach in [?H9] etc., in [12] the authors studied 
the wellposedness of a semilinear parabolic equation subject to a nonlinear dynamical 
boundary condition with subcritical growth for nonlinearities and then investigated the 
large time behavior of the solution such as blow-up phenomenon. In our present paper, 
we are now interested in the question whether the global solution to a class of semilin- 
ear parabolic equation with subcritical nonlinearity and dynamical boundary condition 
fll.ip - fll.3p will converge to an equilibrium as time goes to infinity. 



Before stating our main results we make some assumptions on nonlinearity /. 
(Fl) f(s) is analytic for s£K. 
(F2) 



\f(s)\<c(l + \s\n 



Vs e M, p e [0, a) 



where 




(F3) 



• for jM — 1, 



lim inf 

\s\— s>+oo 



/(*) 



> 




S 



2 



where A > is the best Sobolev constant in the following imbedding inequality 




Vu\ 2 dx + / u 2 dS > A / u 2 dx; 



• for n = 0, 



lim inf 

\s\— >+oo 



> 0. 



s 



Remark 1.1. Assumption (Fl) is made so that we can use an extended Lojasiewicz- 
Simon inequality to prove our convergence result (see Theorem \1. 6 2\ below). Assumption 
(F2) implies that the nonlinear term has a subcritical growth. Under this assumption, we 
are able to prove the local existence and uniqueness of solutions by adopting the argument 
in fTBj with some modifications. Moreover, due to the subcritical growth, existence of a 
strong solution to stationary problem can be obtained by variational method. Assumption 
(F3) is a condition needed to ensure global existence of solution to our problem fll.ljl - 
fll.3p . For simplicity of exposition the nonlinear term f is assumed to depend only on u. 
However, the results in this paper still remain true for the nonlinear term f(x,u) with 
additional smoothness assumption with respect to x. 

Throughout this paper we simply denote the norm in L 2 (Q) by || • || and equip H 1 (Q) 
with the equivalent norm 



We are now in a position to state our main results. 

Theorem 1.1. Suppose (F1)(F2)(F3) are satisfied. Then for any initial data Uq G 
problem 111. i)) -/ TO|) admits a unique global solution such that 




u e c([o, +00); H l {tt)) n c((o, +00); H 2 (tt)) n ^((o, +00); L 2 (n)) 



j(u) e c([o, +00); ffi(r)) n c((o, +00); fft(r)) n ^((0, +00); h?(t)), 



where j(u) is the trace operator. Moreover, if Uq G H 2 



(fl), then 



ueC{[0,+oc);H 2 {{})), 7 (u) eC([0,+oo);#§(r)). 
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Theorem 1.2. Suppose that (Fl) (F2) (F3) are satisfied. Then for any initial data 

4:71 4:71 

Uo G H 1 ^), and in addition uq G L™^(Vt), 7(^0) G L"^(T) when n > 3, the global 
solution to problem U.1\) - [T73\) converges to an equilibrium ij)(x) in the topology of H 1 (Q) 
as time goes to infinity, i.e., 



Here ip(x) is an equilibrium to problem / ti.i| )-/ TOj) . i.e., ip(x) is a strong solution to the 
following nonlinear elliptic boundary value problem: 



Before giving the detailed proof of our convergence theorem, we briefly recall some 
related results in the literature. For the semilinear parabolic equation with homogeneous 
Dirichlet boundary condition, concerning convergence to an equilibrium as time goes to 
infinity, we notice that in one spacial dimension, the situation is much simpler: the set 
of stationary states is discrete and one can take advantage of existence of a Lyapunov 
functional to prove that any bounded global solution will converge to an equilibrium (see 
e.g. |22|I33|). In the case of higher spatial dimension, it was proved in [T7l|26] by use of 
the so-called Lojasiewicz- Simon inequality that any bounded global solution will converge 
to a single stationary state provided that / is real analytic. This kind of result is highly 
nontrivial, since in higher spatial dimension the stationary states can form a continuum 
(see, e.g., [13]). Moreover, a counterexample for semilinear parabolic equations with 
C°° nonlinearities are given in the literature (see [24J, also [23]) to show that there is a 
bounded global solution whose w-limit set is diffeomorphic to the unit circle S 1 . After 
this breakthrough, many further contributions were made for various types of nonlinear 
evolution equations (see, e.g., |2llT0 | [TT | [T4 l [T6 | fT8 | l23 |l29l - [32] and references therein). 

Some new features of our present work are as follows. First, due to the presence of 
dynamical term u t , it turns out that for the corresponding elliptic operator, the dissi- 
pative boundary condition should be considered as non- homogeneous. As a result, the 
tojasiewicz-Simon type inequality we are going to derive is naturally different from the 
one with the usual homogeneous boundary conditions in the literature. Furthermore, 
we are able to treat both non-homogenous Neumann and Robin boundary conditions 
(corresponding to \x = 0, 1, respectively). Second, the subcritical growth of nonlinearity 
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seems natural as stated in Remark 1.1. Besides, we do not have any restriction on spa- 
tial dimension. Thus, our result improves the growth condition of nonlinear term stated 
in jT6l Chapter 3, Section 3.2] even for the Dirichlet boundary condition. 

The remaining part of this paper is organized as follows. In Section 2 we prove the 
existence and uniqueness of the global solution (Theorem II. ip and derive some uniform a 
■priori estimates. In Section 3 we first study the stationary problem. Then we proceed to 
establish the generalized Lojasiewicz-Simon type inequality, and to complete the proof of 
Theorem 11.21 



2 Global Existence and Uniqueness 

Let 

E := {(u,v) T G H 2 (Vl) x H%(T) : j(u) = v}; 

E 1 := {(u,v) T e H\n) x H^(T) : ^(u) = v}; 

F := L 2 (tt) x H*(T). 

Following [7ti9|fT2] we can rewrite problem ( ll.ip - (jl.3p in the following form: 



d_ 

at 



u \ ( u \ ( f(u) 



VI \ V 



) (2.1) 



(«(0)^(0)) t =Kt(«o)) T 



with 

.4 







\ v I 





-An 



and E is the domain of A, i.e., D(A) = E. 

By the results in Escher [9] and Fila & Quittner [12] . we have 

Theorem 2.1. Let uq G iJ 1 (fi) and assumptions (F1)(F2) be satisfied. Then problem 
( fi.-fp -j TOP admits a unique maximal solution u(t) such that 

Z := (it, 7 H) T G C([0, T max )- E t ) n C((0, T max )- E) n ^((0, T max )- F). 

Moreover, the solution u(t) exists globally if u(t) is bounded in H 1 ^). In addition, if 
u G H 2 (n), then Z G C([0,T max ); E) D C\[0,T max ); F) . 



Proof. Let u$ G H l (Q). By assumption (F2), the condition (W) in Escher [9] is satisfied. 
Thus Z G C([0,T max ); Ei) follows from [9j Theorem 1]. On the other hand, it follows 
from [HI Theorem 2.2] that Z G C((0, T max ); £7) n ^((O, T maa; ); F) and furthermore Z G 
C([0, T max ); E) n CHfO, T ma:c ); F), if Mo e # 2 (^). □ 



Define 



E{u{t)) = \ [ \Vu\ 2 dx + ^ I u 2 dS+ I F{u{t))dx (2.2) 



2 Jn 2 j r j n 



where F(z) = f* f(s)ds. 

It's easy to see that V t G (0, T maa .), 

^F(n(t)) + ||ni|| 2 + ||nj2 2(r) = , (2.3) 

which implies that -E(w) is decreasing respect to time. 

Next we prove the global existence of the solution to problem (jl.ip — (jl.3p . 

Theorem 2.2. Assume (F1)(F2)(F3) hold. Then the solution u(t) obtained in The- 
orem 2.1 exists globally, i.e., T max = +oo. Furthermore, assuming in addition that 
Uq G L~^{VL), 7(110) G L^(T) when n > 3, then for any o > we have the follow- 
ing uniform estimate 

\\u\\trz(p) <Ca, V t > a, 

where C a is a positive constant depending only on ||iio||.ffi(n)j ll 7/ oil^ w 4 " 2 ^ . ||7(mq) ||^ , 
a and \Q\. 

Proof. By Theorem 2.1, in order to get the global existence, it suffices to obtain uniform 
a priori estimates on 

From the Sobolev imbedding theorem and the growth assumption (F2), we can get 

F(u)dx = [ [ f(s)ds 
<) Jn Jo 

\u 



< cl I (1 + \s\ p+1 )dsdx < C I {\u\ + \u\ p+2 )dx 
Jn Jo Jn 



< Cfdlulln^ + IHlSitn^CdHl^n)). (2.4) 

This implies 

E(u(t)) < C(\\u\\ HHn) ). (2.5) 

In the same way as in [5], we can also get an estimate in the opposite direction. 
To see this, for sufficiently small 5 > 0, we consider two cases. 
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Case 1: fj, — 1. 

From the assumption (F3), there exists iV = N(5) being a positive constant such that 
f(z)/z > -A + 25 for \z\ > N. Then we have 

Jo z Jn z 2 

for \s\ > (|(iV 2 - 2C)) 1/2 := M, where C = J Q N ^zdz. 

For negative s one can repeat the same computation with N replaced by — N. Now we 
have 

f F{u)dx= [ F(u)dx+ [ F{u)dx> [ u 2 dx + C(\tt\,f) (2.7) 

Jn J\u\<M J\u\>M 2 J n 

where C(|Q|, /) = |Q|minF(s). 

\s\<M 

Thus we can deduce that 

E(u(t)) > e Q j[ \Vu\ 2 dx + X - jf u 2 dS^j + C(\Q\, f) (2.8) 
provided e < ^. Taking C = 1/e, we get 

\\\u\\ 2 m { n) = \ I \Vu\ 2 dx+ l - j u 2 dS<C(E(u(t))-C(\Q\,f)). (2.9) 
Case 2: // = 0. 

By assumption (F3), there exists N = N(5) being a positive constant such that f(z)/z > 
25 for \z\ > N. Then we have 



F(s) = I J -^zdz + / J -^-zdz > -s 2 (2.10) 



for s 2 > i;{N 2 — 2(7). Similar to Case 1, we have 



/ F(u)dx= [ F(u)dx+ [ F(u)dx> 5 - [ u 2 dx + C(\n\J) (2.11) 

in J\u\<M J\u\>M 2 



>|<M J|n|>M 
r 

where C(|Q|, /) = minF(s) - -M 2 \tt\. 

\s\<M 2 

Thus we can deduce that 



E(u(t)) > Q jT |v«| 2 ^ + ^ jT + c(|n|, /) 



(2.12) 



Taking C = 1/5, we get 



1 

2 Ir^iiH 1 ^) 



-;: u|lm^ = i / \Vu\ 2 dx + ]- I u 2 dx < C(E(u(t)) - C(\Q\, /))■ (2-13) 



For both cases, we have 

hWhm < c(E( u (t))-c(\n\j))<c(E( Uo )-c(\n\j)) 

< C(\\u \\ m{n) )-C(\n\,f) 

This implies the following uniform estimate: 

\\u\\ H i m <M, Wte(0,T max ) (2.14) 

where M does not depend on T max . According to Theorem 2.1, this leads to T max = +00, 
i.e., the solution exists globally. 

In order to get the uniform bound of u(t) in H 2 (Q) norm, we are going to get some 
higher-order estimates via a formal argument. However, this procedure can be made 
rigorously within an appropriate regularization scheme (e.g. [35], Chapter 6). 

In what follows we only consider the case n > 3 since the cases n = 1, 2 we can get 
the same result with a simpler proof (see Remark 2.1). 

Again we consider two cases. 
Case 1 (i = 1. Multiplying (II. ip by \u\ p ~~ 2 u, integrating over Q, we get 

\jt\j \ u \ Pdx + j W\ p ds)+{p-l) j \u\ p - 2 \Vu\ 2 dx + j f(u)u\u\ p - 2 dx + j \u\ p dS = 

(2.15) 

Taking p = from assumption (F3), we can deduce that for 5 > sufficiently small 
there exists N = N(5) > such that 

^>(- 4{p - 1] X + 5x) forH>iV, n>3, 
u \ p J 

and result 

f(u)u > (^- 4 ^I - A + SX^j u 2 {ot\u\>N. (2.16) 

(I2.16P implies that 

/ f(u)u\u\ p - 2 dx > C f + ( - 4 ^7 ^ A + ^A ) / \u\ p dx (2.17) 
Jn V P J J\u\>n 



where Cf = Ifilmin f(s)s|s| p 2 . Then we get 

J 1 '[ s |<jv j w 1 1 ta 



~( ! \u\ p dx+ [ \u\ p ds] + 4(P - X) f \V\u\h 2 dx+ [ \u\ p dS 
Pdt \Jn Jr J P 2 Jn Jr 

< ( 4 ^7 - A-5A ) / |m| p c/x + C. (2.18) 



P J Jn 



Using assumption (F3) and applying the imbedding inequality 

/ \Vw\ 2 dx + / w 2 dS > X w 2 dx 
Jn Jt Jn 

to \u\ p / 2 , we get 

j t (j \u\ p dx + J \u\ p dS^j + d (^J \u\ p dx + J \u\ p dS^j <C 2 



Integrating with respect to t yields that 



/ \u\ p dx+ [ \u\ p dS<e~ Clt ( [ 
Jn Jt \Jn 



\u \ p dx+ [ \u \ p ds)+% 
n Jt J w 



Thus, we get 
which implies 



\u\ n - 2 dx < C 



\\f'{u)\\ 



L n (n) 



< C. 



Here C is a constant depending on 
Differentiating (II. ip with respect to t, we get 



in , ||7(«n)ll -4« and \Q\. 



u tt - Au t + f{u)u t = 0. 



Multiplying (12.231) by Ut and integrating over fi, we obtain 
1 d 



Ut II 2 + II u t \\ 2 T 2 (T ,)+ || Vu t || 2 + 112 



+ / f'{u)u 2 dx = 0. 
Jn 



2dt K " ~" ~ l ilL 2 (r)J^ II II -t-||«t|lL2(r) 
It follows from Holder's inequality and the Gagliado-Nirenberg inequality that 



f'(u)u 2 dx 



<\\f'(u)\\ L «\\u t f^ 



L »- 



< G||/ / (Tl)|| L n(||V« t ||||« t || 

< ±\\Vu t \\ 2 + C\\u t \\ 2 . 



Hence, we have 



< 



pt 

t (ll u t f + || u t ||i 2(r) ) + J r (|| Vu* || 2 + || u t ||£ 2(r) ) 
C t T \\u t \\ 2 dT + f (\\u t \\ 2 + \\u t \\h {r) )dr. 



From (12. 3p and (12. 5p we can conclude that for any t > 

/t 
r (|| Vut f + || ut ||l 2(r) ) dr < C(l + t). (2.27) 

Multiplying (I2.23P by — Au t and integrating over f2, then from the boundary condition 
and Holder's inequality we get 



1 d 
Ydt 



(ll Vu 4 f + || u t ||| 2(r) ) + || u tt || 2 2(r) + || Au t 



< || /'(«) II Mt ll L ^ (n) ll Am < II ' ( 2 - 28 ) 

By (I2.22p and the Sobolev imbedding theorem, (I2.28P yields that 

lit (" VMt| ' 2 + H^W)) + H Am *H 2 + II M «H W) ^ 5 II ll' + C II M * II ( 2 - 29 ) 
Multiplying (I2.29P by t 2 and integrating from to t, 

t 2 (ll Vnt f + || u t || 2 2(r) ) +^ r 2 || || 2 rfr + 2 jf r 2 |K|| 2 2(r) dr 

< / r || V« ( || 2 rfr + / r || u t ||| 2(r) dr + 2C [ r 2 \\ u t \\ 2 m{Q) dr. (2.30) 
Jo io Jo 

Combining it with (I2.27P yields that 

t 2 \\u t \\ 2 mm < C(l + t 2 ), Vt>0. (2.31) 

Case 2 /i = 0. 
Similarly we have 



1 d 
p dt 



[ \u\ p dx+ [ \u\ p dS ) + (p - 1) / M p - 2 |Vu| 2 cfe+ / /(M)n|n| p - 2 da; = 0, (2.32) 
in ir / in in 



where p = From (F3), for <5 > small there exists iV = N(5) > such that 

f(u)u > 5u 2 for |w| > N, (2.33) 

which implies 



f(u)u\u\ p - 2 dx >C f + 5 / |w| p c/x > (C f - 8N P \Q\) + 5 / \u\ p dx, (2.34) 
where Cf = |0| min/(s)s|s| p_2 . Then we get 

^ ry n p cfe+ y mWj+c y (jvM2i 2 +| M |p)dx<c. (2.35) 
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This immediately gives (I2.22j) . 

Furthermore we have the following corresponding inequalities to (12.26P (I2.30p : 



t ( || u t || z + || u t \\ z L 2(r) ) + I r || Vm* \\ z dr 



< C I T\\u t \\ 2 dr + J^\\u t \\ 2 +\\u t \\ 2 L2{r) )dT. (2.36) 



t 2 || Vu t || 2 + / t 2 || Au t f dr + 2 / T^lutill^dr 



< 2 / r || Vw 4 f dr + 2C [ r 2 || u t \\ 2 m{n) dr. (2.37) 
Jo Jo 

By using the equivalent H 1 ^) norm of itj, it is easy to see that (I2.3ip holds as well. 
Now for both cases, we consider the following elliptic boundary value problem 

Au = u t + f(u), (23g) 
d v u + fiu + u t | r = 0, 

according to the elliptic regularity theory and the trace theorem we have 

\\u\\ HHn) < C(\\u t + /(«)|| + |k t ||^ (r) + \\u\\) < C(\\u t \\ HH n) + ||/(ti)|| + \W\\), (2.39) 

where C is certain positive constant depending only on Q. 
Combining the fact 

2(n + 2) An 

-A L < for n > 3 

7i-2 n-2 

with the result (I2.2ip . assumption (F2) and Young's inequality yields that 

II /(«) II < C. (2.40) 
By (F2T3T]) . (12~39|) . (12T40|) . for any tr > 0, we obtain 

\\u\\ H 2 m <C a , Vt>a (2.41) 
where C a is a constant depending on Hwollff^Q), HwoU^j^ , llT^o)!!^^^, °" an d 

□ 



Remark 2.1. VFe consider the case n> 3 in the proof when we want to get the estimate of 
\\f{u)\\ and also \\f'(u) \\L n (n)- For n = 1,2, the proof is simpler because we can simply use 
ll^oll-ff 1 ^) t° 9 e t the desired estimates due to the Sobolev imbedding inequality for n = 1,2. 
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3 Proof of Theorem 



T72 



The proof of Theorem 11.21 consists of several steps. 
Step 1. 



From f)2.3p it is clear that the energy E(u(t)) is decreasing with respect to time. On 
the other hand, (12.81) or (I2.12p implies that E(u(t)) is bounded from below. Thus we 
know that E(u(t)) serves as a Lyapunov functional. Besides, V t > , if E(u) = E(uq) 
then t| = 0. This enables us to deduce that Ut = 0, which means u is an equilibrium. 
The a»-limit set of Uq G H 1 (Q) is defined as follows: 

u(u ) = {ij(x) eH\Q) : 3 {t n }™ =1 , t n +oo s. t. u(x,t n ) -> ip(x) in H\Q)}. 

Then from the well-known results in the dynamic system (e.g., [27J, Lemma 1. 1.1) it is 
easy to see that 

Lemma 3.1. The uj-limit set of Uq is a non-empty compact connected subset in H 1 ^). 
Furthermore, (i) it is invariant under the nonlinear semigroup S(t) defined by the solution 
u(x,t), i.e, S(t)u(uo) = u>(uq) for all t > 0. (ii) E(u) is constant on u(uq). Moreover, 
co(uq) consists of equilibria. 

Step 2. In this step, we collect some results on the stationary problem. 

Lemma 3.2. Suppose that ip G H 2 (Q) is a strong solution to problem M.5\) . Then ip is a 
critical point of the functional E(u) in H 1 ^). Conversely, if if} is a critical point of the 
functional E(u) in H 1 ^), then t/j G H 2 (n), and it is a strong solution to problem U.5\) . 



Proof. If if) G H 2 {VL) satisfies ( TO]) , then for any v G H 1 ^), it follows from (TT3|) that 



Thus, ifi is a critical point of E(u). Conversely, if if) is a critical point of E{u) in H l (VL), 
then (13. ip is satisfied. By the assumption (F2) on subcritical growth and the bootstrap 




By integration by parts and the boundary condition in (ll.5p . we get 




(3.1) 



which, by a straightforward calculation, is just the following 

dE(ifj + ev) 



argument, if) G H 2 



(Q) follows. 



□ 
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The following lemma claims that problem 01 .5p admits at least a strong solution. 
Lemma 3.3. The functional E(u) has at least a minimizer v G H 1 ^) such that 

E(v) = inf E(u). 
In other words, problem U.5\) admits at least a strong solution. 

Proof. From Section 2 we can see that E(u) is bounded from below on Therefore, 
there is a minimizing sequence u n G H l (Q) such that 

E(u n ) -> inf E(u). (3.2) 

ueH 1 



On the other hand, E(u) can be written in the form: 

E(u) = ^\\Vu\\ 2 + ^\\u\\l 2{r) +F(u) (3.3) 

with 

JF( U ) = / F(u)dx, (3.4) 
Jn 

it follows from (12. 8p ( (I2.12p for [i = 0) that u n is bounded in H l (Q). It turns out from the 
weak compactness that there is a subsequence, still denoted by u n , such that u n weakly 
converges to v in if x (f2). Thus, v G if 1 (f2). We infer from the Sobolev imbedding theorem 
that the imbedding H 1 ^) ^ L 7 (f2) ( 1 < 7 < is compact. As a result, u ra strongly 
converges to v in L 7 (f2). It turns out from the assumption (F2) that J"(u„) — > ^(v). 
Since is weakly lower semi-continuous, it follows from ( 13. 2 p that E(v) = inf E{u). 

From the assumption (F2) on subcritical growth, by the bootstrap argument and elliptic 
regularity theorem, this weak solution is also a strong solution in H 2 (Q). Furthermore, we 
can also use the bootstrap argument to show that v G L°°(Q). The proof is completed. □ 

Remark 3.1. Under assumption (F2), by a further bootstrap argument, we can show that 
ij) is also a classical solution. 

Step 3. Define 

V := {u G H 2 {Q) I d v u + fiu | r = 0}. 

Let ^ be a fixed critical point of E(u). In what follows we are going to establish 
the generalized Lojasiewicz-Simon inequality which extends the original one by Simom 
[26J for the second order nonlinear parabolic equation subject to the Dirichlet boundary 
condition. Similar inequalities have been derived in our recent papers [30-32J to prove 
the convergence to equilibrium for various evolution equations (systems) subject to the 
dynamical boundary conditions. 
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Lemma 3.4. Let ip be a critical point of E(u). Then there exist constants 9 G (0, |) and 
(3 > depending on if) such that for any u G H 2 (Q) satisfying \\u — ipW^tn) < P and 
\\d u u + hu\\l2( T } < [3, we have 

\\-Au + f(u)\\ (H i m + \\d v u + ^u\\ L 2 {r) > \E(u) - E{^)\ 1 - 6 . (3.5) 

Proof. Let 

M(u) = -Au + f(u). (3.6) 
Then M(u) maps u G H 2 (Q) into L 2 (Vt). 

First, we claim that there exist constants 9' G (0, |) and 0i > depending on if), Vu; G 

T> (i.e., satisfies the homogeneous Neumann or Robin boundary condition) satisfying 
|| w — ip\\m < Pi, we have 

\\M(w)\\ {HHn)Y > C\E(w) - Em 1 '" ■ (3-7) 

The above claim follows from the same argument as that in |14[[TT] for the case of homo- 
geneous Dirichlet boundary condition (see also |2S]), so the details can be omitted here. 
However, we are now dealing with the nonhomogeneous boundary condition. As a result, 
we have to extend the above result from functions in T> to H 2 (Q). Hereafter we always 
keep in mind that if) is an equilibrium satisfying d v ip + flip = on T. As before we consider 
two cases. 

Case 1 fx = 1. For any u G H 2 (Q), we consider the following Robin Problem: 

Aw = Au, x G Q, 
d u w + w = 0, x G r. 

It has a unique solution w belonging to D. On the other hand, we introduce the Robin 
map R : H S (T) — > H s+ ( 3 / 2 \tt) which is defined as follows: 

f Aq = in Q, 
Rp = q < 

[ d u q + q = p on T. 

As shown in [H], R is continuous for s G M. Thus we have 

\\w - u\\ H i(Q) < C\\d u u + wilier), (3.8) 
where the constant C does not depend on u. As a result, 

\\w - ^lUi(n) < \\w - u\\ H i(n) + ||w - 4>\\m(n) < C\\d u u + u\\ L 2 {v) + ||w - ip\\w-(ay (3-9) 
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Then it is easy to see the fact that u G H 2 (Q) being in the small neighbourhood of ip in 
and \\d u u + u\\l2^ being small imply that w stays in a small neighbourhood of ip 
in H\Sl). 

By direct calculation and the Sobolev imbedding theorem we can see that 



\\M(w)\\ imm , < \\M{u) 

< \\M(u) 

< \\M(u) 

< \\M(u) 

< \\M(u) 



(HHn))> + \\f(u) - f(w)\\ iH i m > 



f'{u + t{w — u)){u — w)dt 



(HHQ))' + max \\f(u + t(w - u))\\ l9 {a) \\u - w\\ L ^ {n) 
(Hi(Q))' + C\\u — w\\h^{0.) 

(fli(n))' + C\\d v u + u\\ L 2 {T) . (3.10) 



On the other hand, by the Newton-Leibniz formula and (13.81) we can get 



\E(w) - E(u) 



< 



JO, 



M{u + t{w — u)){u — w)dxdt 



+ 



o Jr 



1 — t)(d u u + u){u — w)dSdt 



< C (||M(«)||(£ri(n))' + \\duU + u\\ L 2 {T) ) \\d u u + u|| L 2 (r) 

< C (||M(«)||(£ri(n))' + \\dvu + M||L2 (r) ) 2 . 



Since 



E(w) - E{$) \ x - e '>\ E(u) - E(ip) 



1 1—9' 



E(w)-E(u) \ 1 - e \ 



(3.11) 



(3.12) 



and < & < §, 2(1 - 6') - 1 > 0, then from ([SID and (l3~T0|) - (l3~TT] ) we can conclude that 

C(||M(n)|| (// i ( n)y + ||^n + n|| L2(r) ) > |£(u) - E^ 1 ' ', 

for 1 1 it — ^llfmfi) < (3, \\d u u + m|| L 2(p) < (3, where (3 > is chosen small enough such that 
\\w - i/j\\ h i < ft. 

Next we choose e, < e < 6' and (3 smaller if necessary such that when ||w — ipWn 1 < (3, 

1 



C 



E(u) - E(i))\- £ > 1. 



(3.13) 



Setting 6 = 9' — e G (0, |), when ||u — V'llfp < /3 and \\d v u + u\\i^(r) < (3, we finally have 
\\M(u)\\ {H i m , + \\d v u + u\\ LHT) > lE^-E^l 1 - 9 . (3.14) 
which is exactly (13.51) . 

Case 2 /z = 0. For any u G H 2 (Q), we consider the following Neumann Problem: 

—Aw + w = —Au + u, x G Q, 
d u w = 0, x G r. 
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It's easy to see that there exists a unique solution w belonging to T> and 

IMItficn) ^ C (IMIff^n) + IIM| L 2 (r) ) . 

Again we can see that u being in the small neighbourhood of ip in with ||c^tt||z,2( r ) 

small implies that w stays in a small neighbourhood of ip in 
Furthermore, by the energy estimate, we have 



\\w - u\\ H i (n) < C\\d v u\\ L 2 {v) . 
In a way similar to Case 1, the following estimates follow: 



(3.15) 





\\M{w) 


||(.ffi(o))' 




< 


\\M{u) 


l(Hi(n))' 


+ \\(f(u)-u)-(f{w)-w)\\ ma) y 


< 


\\M(u) 


\(HHQ))> 


+ max II f'( u + f ( w ~ U ))W l9( ) II n " 


< 


\\M(u) 




+ C||w - 


< 


\\M(u) 




+ C\\d u u\\ L 2 {T) . 



2n + U 



Hlc^n))' 
(3.16) 



On the other hand, by the Newton-Leibniz formula and (I3.15p . we can get 



< 



E(w) - E{u)\ 

I M(u + t(w — u))(u — w)dxdt 
o Jn 



+ 



o Jr 



1 — t)d u u{u — w)dSdt 



< C(||M(u)|| (HW + ||^|| i2(r) ) . 



(3.17) 



Then arguing as in Case 1, we can get the required corresponding result for fi = 0. 

Thus, the lemma is proved. □ 

Following the idea in |32] (see also j3U]) we can easily extend the previous lemma in 
such a way that we only need the smallness of \\u — ipW^tn), 

Lemma 3.5 (Generalized Lojasiewicz-Simon Inequality). Let ip be a critical point of 
E(u). Then there exist constants 9 G (0,^) and /3 > depending on ip such that for 
Vu G H 2 {Vt), \\u — i>\\m(n) < Ab we have 



- Au + f{u)\\ {mm , + \\d u u + ^u\\ L 2 {r) > \E(u) - E(il))\ l ~ e . 



(3.18) 
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Proof. Let (3 > and 9 G (0, |) be the constants appearing in Lemma [3.41 
We consider the following two cases: 

(i) ||w — ipWn 1 ^) < ft and \\d u u + /i-uH^p) < (3. Then from Lemma [3.41 we are done. 

(ii) \\u — 1^11^(0) < P but \\d u u + fm\\ L 2^ > (3. By Sobolev imbedding e — >■ L^(f2), 
there exist /3 > depending on t/> such that for any u satisfying \\u — V-'IIh 1 ^) < /3, 



-E^l 1 - 6 < (3. (3.19) 



Thus, we have 



|| - Au + f(u)\\ {H i {n)y + \\d u u + fiu\\ L 2 {r) > \\d u u + flu\\ L 2 {r) 

> a 

> \E(u) - E(i))\ 1 ' 6 . (3.20) 

Finally by setting /3 = min{/3,/3}, the lemma is proved. □ 

Step 4. After the previous preparations, we now proceed to finish the proof of Theorem 
II. 2\ following a simplified argument introduced in |17| in which the key observation is 
that after certain time to, the solution u will fall into the small neighborhood of a certain 
equilibrium ip{x), and stay there forever. 
First, it is easy to see that 

||ttt|| —7-0 as t — > +oo. 
From Step 1, there is a sequence {t n } ne N, t n +oo such that 

u{x,t n ) ip{x), mH l (Q), (3.21) 

where ip(x) G u>(uq) is a certain equilibrium. On the other hand, it follows from (12. 3ft that 
E(u) is decreasing in time and 

lim E(t n ) = E(ip). (3.22) 



n— ¥+oo 



We now consider all possibilities. 

(1) . If there is a to > such that at this time E(u) = E(ip), then for all t > t , we 
deduce from (12. 3p that u is independent of t. Since u(x,t n ) — > ip, the theorem is proved. 

(2) . If there is to > such that for all t > to, u satisfies the condition of Lemma [3.5[ 
i.e., \\u — iP\\h 1 (q,) < 00, then for 9 G (0, \) introduced in Lemma l3~5l we have 

±(E(u) - Em 6 = 0(E(u) - EWr 1 ^- (3.23) 
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Combining it with (1515)1 . Q yields 

j t {E(u) - E(if>)) e + \ (II -A« + /(«) || + || |U 2(r) ) < 0. (3.24) 
Integrating from t to t, we get 

(E(u) - E^)) e + lf (|| —Aw + /(«) || + || || L2(r) ) < (E(u(t )) - E(^)) e . (3.25) 
Since -E(m) — ^(VO > 0, we have 

/ (|| -Au + f(u) || + || u t \\ L z(r))dT < +oo, (3.26) 

J to 

which implies that for all t > to, 

ft 

\u t \\dr < C. (3.27) 



to 

This easily yields that as t — > +oo, u(x,t) converges in L 2 (Q). Since the orbit is compact 
in H l (Q), we can deduce from uniqueness of limit that (11.41) holds, and the theorem is 
proved. 

(3). It follows from (I3.2ip that for any e > with e < (3q, there exists an integer N 
such that when n > N, 

II U(;t n ) - V ||< \H;t n ) - lp\\ H i m < |, (3.28) 

l -{E{u{t n ))-Em e <- A . (3.29) 

Define 

t n = sup{ t > t n | || u(-,s) -ip \\m(n)< 0o, Vs e (3.30) 

It follows from (I3.28P and continuity of the orbit in H 1 ^) for t > that f n > £ n for all 
n > N. Then there are two possibilities: 

(i) . If there exists no > N such that t no = +oo, then from the previous discussions in 
(1)(2), the theorem is proved. 

(ii) Otherwise, for all n > N, we have t n < i n < +oo, and for all t G [t n ,i n ], E(ip) < 
E(u(t)). Then from (I3.25P with t being replaced by t n , and t being replaced by t n we 
deduce that 

I ut || dr < ^(E(u(t n )) - E{i,)f < |. (3.31) 
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Thus, it follows that 



u{t n ) -ip\\ < \\u(t n ) -lf>\\ + 




u t || dr < e 



(3.32) 



which implies that when n — > +00 



u{t n ) -> if), 



in L 2 (n). 



(3.33) 



Since [J t>s u(t) is relatively compact in H 1 (Q), there exists a subsequence of {u(t n )}, still 
denoted by {u(t n )} converging to if> in H l (Q). Namely, when n is sufficiently large, 



which contradicts the definition of t n that ||w(-,t n ) — V'll^fn) = A)- 

Thus, Theorem 11.21 is proved. 
Acknowledgement: The author is indebted to Prof. S. Zheng for his enthusiastic help 
with this paper. This work is supported by the NSF of China under grant No. 10371022, by 
Chinese Ministry of Education under the grant No. 20050246002, by the Key Laboratory 
of Mathematics for Nonlinear Sciences sponsored by Chinese Ministry of Education and 
the Graduate Student Innovation Foundation of Fudan University. 

References 

[1] R. Adams, Sobolev Spaces, Pure and Applied Mathematics series; a series of mono- 
graphs and textbooks, 65, 1970. 

[2] S. Aizicovici, E. Freireisl and F. Issard-Roch, Long-time convergence of solutions to 
a phase-field system, Math. Methods in the Appl. ScL, 24, (2001), 277-287. 

[3] J.M. Arrieta, P. Quittner and A. Rodriguez-Bernal, Parabolic problems with non- 
linear dynamical boundary conditions and singular initial data, Diff. and Int. Eqns, 
14, No.12 (2001), 1487-1510. 

[4] G.I. Bizhanova and V.A. Solonnikov, Solvability of an initial-boundary value problem 
with time derivative in the boundary condition for a second parabolic equation in a 
weighted Holder function space, St. Petersburg Math. J. 5, No. 1, (1994), 97-124. 

[5] I. Chueshov, M. Eller and I. Lasiecka, On the attractor for a semilinear wave equation 
with critical exponent and nonlinesr boundary dissipation, Comm. Part. Diff. Eqns., 



u(t n ) - ^Hffi(n) < A), 



(3.34) 



27(2002), 1901-1951. 



19 



[6] R. Chill, On the Lojasiewicz-Simon gradient inequality, J. Func. Anal., 201(2003), 
572-601. 

[7] J. Escher, Quasilinear parabolic systems with dynamic boundary conditions, Comm. 
Part. Diff. Eqns. 18(1993), 1309-1364. 

[8] J. Escher, On quasilinear fully parabolic boundary value problems, Diff. and Int. 
Eqns. 7(1994), 1325-1343. 

[9] J. Escher, On the qualitative behavior of some semilinear parabolic problems, Diff. 
and Int. Eqns. 8(1995), 247-267. 

[10] E. Feireisl, F. Issard-Roch, and H. Petzeltova, A non-smooth version of the 
Lojasiewicz-Simon theorem with applications to nonlocal phase-field systems, J. Diff. 
Eqns, 199(2004), 1-21. 

[11] E. Feireisl and F. Simondon, Convergence for semilinear degenerate parabolic equa- 
tions inseveral space dimensions, J. Dyn. Diff. Eqns., 12, No. 3, (2000), 647-673. 

[12] M. Fila and P. Quittner, Large time behavior of solutions of a semilinear parabolic 
equation with a nonlinear dynamical boundary condition. Progress in Nonlinear Diff. 
Equ. and Their Appl. 35, (1999), 251-272. 

[13] A. Haraux, Systemes dynamiques dissipatifs et applications, Masson, Paris, 1991. 

[14] A. Haraux and M.A. Jendoubi, Convergence of bounded weak solutions of the wave 
equation with dissipation and analytic nonlinearity, Calc. Var. 9, (1999), 95-124. 

[15] T. Hintermann, Evolution equations with dynamic boundary conditions, Proc. Roy. 
Soc. Edinburgh, Sect A, 113 (1989), 43-65. 

[16] S-Z. Huang, Gradient Inequalities, with applications to asymptotic behavior and sta- 
bility of gradient-like systems, Mathematical Surveys and Monographs, 126, AMS, 
2006. 

[17] M.A. Jendoubi, A simple unified approach to some convergence theorem of L. Simon, 
J. Func. Anal, 153 (1998), 187-202. 

[18] M.A. Jendoubi, Convergence of global and bounded solutions of the wave equation 
with linear dissipation and analytic nonlinearity, J. Diff. Eqns., 144 (1998), 302-312. 



20 



[19] Joachim von Below and Gaelle Pincet Mailly, Blow up for reaction diffusion equations 
under dynamical boundary conditions. Comm. Part. Diff. Eqns., 28 (2003), no. 1-2, 
223-247. 

[20] W. Li, Long-time convergence of solution to phase-field system with Neumann bound- 
ary conditions, to appear in Chinese Ann. Math. A. 

[21] S. Lojasiewicz, Une propriete topologique des sous-ensembles analytiques rees. Col- 
loque Internationaux du C.N.R.S. 117, Les equations aux derivees parielles (1963), 
87-89. 

[22] H. Matano, Convergence of solutions of one-dimensional semilinear parabolic equa- 
tions, J. Math. Kyoto Univ. 18-2 (1978), 221-227. 

[23] P. Polacik and K.P. Pybakowski, Nonconvergent bounded trajectories in semilinear 
heat equations, J. Diff. Eqns., 124 (1996), 472-494. 

[24] P. Polacik and F. Simondon, Nonconvergent bounded solutions of semilinear heat 
equations on arbitrary domains, J. Diff. Eqns., 186 (2002), 586-610. 

[25] P. Rybka and K.H. Hoffmann, Convergence of solutions to Cahn-Hilliard equation, 
Comm. Part. Diff. Eqns., 24 (5&6), (1999), 1055-1077. 

[26] L. Simon, Asymptotics for a class of nonlinear evolution equation with applicationa 
to geometric problems, Ann. of Math., 118 (1983), 525-571. 

[27] R. Temam, Infinite- dimensional Dynamical Systems in Mechanics and Physics, Appl. 
Math. Sci., 68, Springer- Verlag, New York, 1988. 

[28] H. Triebel, Interpolation Theory, Function Spaces, Differential Opeartors, North- 
Holland, Amsterdam, 1978. 

[29] H. Wu, M. Grasselli and S. Zheng, Convergence to equilibrium for a parabolic- 
hyperbolic phase-field system with Neumann boundary conditions, to appear in 
M3AS, 2006. 

[30] H. Wu, M. Grasselli and S. Zheng, Convergence to equilibrium for a parabolic- 
hyperbolic phase-field system with dynamical boundary condition, to appear in J. 
Math. Anal. Appl., 2006. 

[31] H. Wu and S. Zheng, Convergence to equibrium for the Cahn-Hilliard equation with 
dynamic boundary conditions, J. Diff. Eqns., 204, (2004), 511-531. 



21 



[32] H. Wu and S. Zheng, Convergence to equilibrium for the damped semilinear wave 
equation with critical exponent and dissipative boundary condition, Quart. Appl. 
Math. 64 (2006), no. 1, 167-188. 

[33] T.I. Zelenyak, Stabilization of solutions of boundary value problems for a second- 
order parabolic equation with one space variable, Differentsial'nye Uravneniya, 1968, 
17-22. 

[34] Z. Zhang, Asymptotic behavior of solutions to the phase-field equations with Neu- 
mann boundary conditions, 4 (2005), no. 3, 683-693. 

[35] S. Zheng, Nonlinear Evolution Equations, Pitman Monographs and Surveys in Pure 
and Applied Mathematics, 133, CHAPMAN & HALL/CRC, Boca Raton, Florida, 
2004. 



22 



